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The dynamics of fluid vesicles in simple shear flow is studied using mesoscale simulations of 
dynamically-triangulated surfaces, as well as a theoretical approach based on two variables, a shape 
parameter and the inclination angle, which has no adjustable parameters. We show that between 
the well-known tank-treading and tumbling states, a new "swinging" state can appear. We predict 
the dynamic phase diagram as a function of the shear rate, the viscosities of the membrane and the 
internal fluid, and the reduced vesicle volume. Our results agree well with recent experiments. 
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The dynamics of deformable objects such as liquid 
droplets lEJl, lipid vesicles [3, 0, i, d, 0, i, 0, [H, red 



13,ll3|,llJ]j and elastic capsules [l 




blood cells [11 

in flows has received increasing attention experimentally, 
theoretically, and numerically in recent years. All of 
these objects show phenomenologically similar behaviors 
in shear flows — either a tank-treading rotation with 
a stationary shape and a finite inclination angle 9 > Q 
with respect to the flow direction, or an unsteady tum- 
bling motion. However, the qualitative and quantita- 
tive behavior of the various objects can be very different, 
because the energies governing their shapes and ther- 
mal fluctuations are very different. In the case of fluid 
vesicles, which we want to investigate in this letter, the 
membrane conformations are determined by the curva- 
ture elasticity together with the constraints of membrane 
incompressibility and constant internal volume. 

In simple shear flows, with flow velocity v = jye^, a 
transition occurs from tank-treading to tumblingwith an 
increasing viscosity of the internal fluid Tyin 0, M or 
membrane viscosity rjmh @, This transition is de- 

scribed very well by the theory of Keller and Skalak (K- 
S) [l3|, which assumes a flxed ellipsoidal vesicle shape. 
However, shear forces can be large enough to induce 
shape transformations of fluid vesicles, for example from 
discocyte to prolate at small membrane viscosity and 
small viscosity contrast between inside and outside, or 
fro m p rolate to discocyte at larger membrane viscosities 
UUM- In this case, it is essential to take the membrane 
deformability into account. 

Recently, Kantsler and Steinberg [3| reported the flrst 
observation of a new type of vesicle dynamics in shear 
flow, which is characterized by small oscillations of the 
inclination angle 9 and the deformation, where —9q < 
(^{t) ^ ^0 with 9o/n < 1 and time average {9) ~ 0. The 
vesicles were found to transit from tumbling to this os- 
cillatory motion with increasing shear rate 7. It is worth 
mentioning that such an oscillatory motion was also ob- 
served in our previous simulations (see Fig. 6 in Ref . [loh . 
but not further analyzed. Simultaneously with the ex- 
periment, Misbah Q predicted a "vacillating breathing" 
mode for quasi-spherical fluid vesicles. This mode ex- 



hibits very similar dynamical behavior as seen experi- 
mentally; however, it "coexists" with the tumbling mode, 
and its orbit depends on the initial deformation i.e. it is 
not a limit cycle. Furthermore, the shear rate only enters 
the theory as the basic time scale, and therefore cannot 
induce any shape transitions. Hence it does not explain 
the tumbling-to-oscillatory transition seen in the experi- 
ments |4|. 

In this letter, we study the oscillatory dynamics of fluid 
vesicles, which we will refer to as the "swinging mode" , 
by mesoscale hydrodynamics simulations and a simpli- 
fled non-linear theoretical model. The main questions 
we want to address are: How does the bending energy 
affect the dynamics? Can transitions between modes be 
induced by varying the shear rate? What happens be- 
yond the quasi-spherical limit, which is the typical exper- 
imental situation? What is the effect of the membrane 
viscosity? Is the swinging mode stable when thermal 
membrane undulations are taken into account? We will 
show that the experiments of Ref. ^ can be understood 
very well on the basis of our theory. Furthermore, we will 
present a complete phase diagram of vesicle motion as a 
function of shear rate and viscosity. 

It is worth mentioning that elastic capsules [l5| and 
red blood cells [l3| can also exhibit a swinging mode; 
however, in this case, the angle 9{t) is always positive 
during the oscillation. Very recently, this dynamics was 
explained by the K-S theory with an addition of an en- 
ergy barrier for the tank-treading rotation caused by the 
membrane shear elasticity Therefore, this mecha- 

nism cannot be employed to explain the swinging mode 
of fluid vesicles. 

The vesicle dynamics is described by several dimen- 
sionless quantities. For a vesicle of volume V and sur- 
face area 5", the reduced volume V* and the excess area 
As are deflned by V* = {Ry/Rsf = (1 + As/An)-^/^ 
and As = S/Rl - Att, where i?v = {W/At:)^^^ and 
Rs — {S/Any/^. The relative viscosity of the inside fluid 
and membrane are 77^^^ = Vin/m and 77^"^^ = ?7mb/?7o-Rs, 
where 770 is the viscosity of the outside fluid. The shape 
relaxation time of vesicles with bending rigidity n is given 
by T = r]oRg/K (for 77*^^ = 1). This time is used to deflne 
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FIG. 1: (Color online) Temporal evolution of au and inclina- 
tion angle 8, for V* = 0.78 and Ty^b ~ 2.9. The solid (red) and 
dashed (blue) lines represent simulation data for 7* = 3.68 
and 0.92 {K/ksT = 10 and 40 with -yrioRl/kuT = 36.8), re- 
spectively. The solid lines in (b) are obtained from Eqs. ((6}, 
(O without thermal noise for 7* = 1.8, 3 and 10 (from top to 
bottom). 



a reduced shear rate 7* — 7r. 

The hydrodynamics of fluid vesicle can be studied very 
well by hybrid simulations of a dynamically-triangulated 
membrane model J^j and a particle-based mesoscale sol- 
vent, multi-particle collision dynamics (20l |. A detailed 
description of this method can be found in Ref. [l^. 

We have simulated a prolate vesicle with a viscous 
membrane at V* = 0.78 and rj^^ — 1. Figure [1] shows 
the time development of the vesicle shape and 9. The 
shape parameter is a-Q = (Li — L2)/(ii + L2), where Li 
and L2 are the maximum lengths in the direction of the 
eigenvectors of the gyration tensor in the vorticity {x^y) 
plane. The vesicle is found to exhibit a tumbling motion 
for 7* — 0.92, but a swinging motion for 7* = 3.68. As 
7* increases, the tumbling frequency /tmb decreases, see 
Fig. [21 while the frequency of the swinging motion in- 
creases. The shear rate 7* at the tumbling-to-swinging 
transition is found to increase with increasing membrane 
viscosity These results are consistent with the ex- 

periments of Ref. 0- 

In order to obtain a better understanding of the vesi- 
cle dynamics, we now derive an ajpproximate theoreti- 
cal model. First, we follow Refs. 0, and employ the 
Stokes approximation and perturbation theory for quasi- 
spherical vesicles. The vesicle shape is expanded in spher- 
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FIG. 2; (Color online) 7 dependence of the tum- 
bling frequency /tum normalized by the frequency /t^m = 
7 VI - of K-S theory at V* = 0.78. The bending rigid- 

ity K is varied to change 7* with •yrfoR^/k-BT — 36.8. The 
symbols represent simulation data for r]^^ = 2.9 (•) and 4.1 
(■). The solid and dashed lines are obtained from Eqs. ([Sjl, 
with and without thermal noise, respectively. 



ical harmonics Yi^rn as i? = Ry{l + J2i,m ui.mXum)- The 
free energy of the membrane with bending rigidity k and 
surface tension cr is F = / dS{(j + (k/2)(Ci -I- C2Y}, 
where Ci and C2 are the principal curvatures at each 
point of the membrane. With the harmonic approxima- 
tion, it is given by F = (k/2) Y.l„i Ei\ui^rn? + 0(u/,m^), 
where Ei = [l + 2){l - 1){1{1 + 1) + uRy^/k}. The flow 
fields inside and outside of the vesicle are described by 
the Lamb solution with no-slip boundary conditions 
on the membrane. The flow stress on the membrane is 
balanced with the elastic forces due to bending and ten- 
sion. This implies that the undulation amplitudes 
are determined by 



l.rn 



kT,E, 



dt 



-ui. 



■m , ± 2 ; 



(1) 



where /i = 60^27r/15/(32 + 23?7*„) and F, = + + 
2){2P ^l + 2) + {l- 1){2P + 51 + 5)?7*„}. The detailed 
derivation of Eq. ([1]) for tjI^^ = 1 is described in Ref. 5'. 

In this letter, we focus on vesicles, which are symmetric 
with respect to the vorticity plane through their center. 
Hence, we only consider U2.±2, "2,01 a-nd decompose U2,±2 
into amplitude and phase, U2,±2 = rexp(^2i0), where 9 
corresponds to the inclination angle [G]. Since the cur- 
vature energy does not contribute in this case, compare 
Ref. @, we have to go beyond the harmonic approxima- 
tion. Therefore, we replace the force kEit in Eq. |T]) by 
dF/dr. This implies 
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FIG. 3: (Color online) 7 dependence of the amplitudes of the 
oscillations of qd and 9 for V* = 0.9, rj*^ = 6, and j;^^ = 0- 
The dashed and solid lines represent the tumbling and swing- 
ing states, respectively, given by Eqs. (|6} , (O without thermal 
noise. Symbols represent estimates from the experimental 
data Q for ao (A) and 6 (♦). The inset shows the curva- 
ture energy _F(q:d) of an ellipsoid for the reduced volumes 
V = 0.78,' 0.85, and 0.9. 



A Lagrange multiplier A = (jRy is employed to satisfy 
the area constraint [2l|; it is determined by dAs/dt = 0, 
which implies As = 2u2,o^ -f 4r^. Thus, Eq. ([2]) becomes 



dr ( 4r^ 
'dt^X 



r, dF 



2770 -Ry '^'^ 



A. 



+7/isin(26') L (5) 



where dF/dr\As — dF/dr — {2r/u2fl)dF/du2fl- In the 
harmonic approximation of F , the first term in Eq. 
disappears because of dF/dr\As = 0, so that we recover 
the description of Ref. @. The prolate shape appears as 
an energy minimum, when higher-order terms in the free 
energy F are taken into account 21]. 

Since an expansion of the vesicle shape in spherical 
harmonics is difficult experimentally, more easily acces- 
sible measures of the deformation are desirable. Fur- 
thermore, a description is needed which goes beyond 
the quasi-spherical limit. Thus, instead of r, we em- 
ploy the shape parameter — (Li — L2)/{Li + L2), 
which is easily measurable by microscopy. Since — 
(^15/27r/2)r -KO(r2), p implies 




Aq dF 
tkV* dau 



jAi sm{29) 



(6) 

where ^0 = ISFa/lGTr = 45/87r(32 + 23r]*J and Ai = 
/iA/l5/27r/2 = 30/(32 -I- 23r/*J. Since an accurate eval- 
uation of the free energy F is very important, we calcu- 
late it numerically for ellipsoidal vesicles with (xi/ai)^ + 
{x2/a2Y + [x^/a^Y = see the inset in Fig. O The 
prolate (oi > 02 = a^) and oblate (oi = a2 > 03) 
shapes are energy minima and maxima, respectively, and 
dF/daD diverges in the limit of maximum extension, 
ao af^^^iV*). Eq. ((6|) has the same form as the sim- 
plified model studied previously [^, [lo| , but now has no 
adjustable parameters. 



In a final step, in order to obtain a reliable descrip- 
tion also for large excess areas, we replace Eq. ([3]) by the 
equation of K-S theory I3, [3 j which reads 



^ - i{-l + i^cos(2.)}. 



(7) 



B - /o ./i + 



l + /2(?7rn-l) + /2/3<bJ ' 



where the factors /o, /i, /2, and /a are the functions 
of the ellipsoidal shape (02/01 , 03/01), and are given 
in Appendix B of Ref. |lO|. The K-S theory in general 
shows very good agreement with simulation results 0, 
[lo| . When thermal fluctuations are taken into account, 
Gaussian white noises ga{t) and geit) are added to 
Eqs. ([6]) and ([7]), respectively, which obey the fluctuation- 
dissipation theorem, so that {gi{t)) = 0, {gi{t)gj{t')) = 
2kBT/(idijS{t — t'), where i,j G {au,d}. The friction co- 
efficients are = tkV* / Ao{l - {a^ / a'^'^'^)'^} and the ro- 
tational friction coefficient of a sphere, C,g — SnrjoRg. We 
numerically integrate Eqs. ((6|) and ([7]) with or without 
thermal noise using the second or fourth-order Runge- 
Kutta method, respectively. When anit) becomes zero, 
we set (aD,^)=(0, 7r/4), since 9 = tt/A is the inclination 
angle in the spherical limit. 

Eqs. ([6]) and ([7]) reproduce the simulation results semi- 
quantitatively (see Figs. [1] and [2]). This is a very good 
agreement, given the fact that our model certainly does 
not systematically take into account all higher order 
terms, and furthermore has no adjustable parameters. 
At small 7*, in the tumbling phase, 9 rotates with au- 
oscillation of small amplitude. This amplitude be- 
comes larger at larger 7*, see Fig. [3l Then, when a^it) 
reaches zero, 9 jumps to 7r/4, see Fig. [H This type of 
9 oscillations with a jump is reminiscent of the behavior 
predicted for very viscous liquid drops within perturba- 
tion theory Finally, at even larger 7*, ao and 9 
exhibit oscillations without jumps. 

The physical mechanism of swinging can be under- 
stood on the basis of Eqs. ©, ([7]) as follows. At finite 
K, the shear force elongates the vesicle {au increases) 
for < < 7r/2, but compresses it (ao decreases) for 
— 7r/2 < 9 < 0, since the sin(26')-term in Eq. ^ changes 
sign. Thus, the swinging motion is caused by a shape de- 
formation, where B in Eq. ([7|) crosses the tank-treading- 
to-tumbling threshold periodically. First, a prolate vesi- 
cle starts tumbling because B < 1, q;d decreases when 
9 < 0, which implies that B increases; then 9 starts to 
increase again because -B > 1 at small q;d; finally q;d 
increases when 9 > 0. In the swinging phase, the ampli- 
tudes of QfD and 9 decrease and saturate to finite values 
with increasing 7*, see Fig. [3] 

A linear stability analysis of the fixed points of Eqs. (O, 
([7|) shows that the tank-treading-to-tumbling transition 
always occurs as a saddle-node bifurcation at 6* ~ 0. The 
tank-treading-to-swinging transition is a saddle-node bi- 
furcation for small 7* , but becomes a Hopf bifurcation at 
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FIG. 4: (Color online) Dynamical phase diagrams as a func- 
tion of (a) for ryf^ = 1, and (b) 77;*^, for r;j^i, = 0, for vari- 
ous reduced volumes V* , obtained from Eqs. ([6l, (O without 
thermal noise. The tank-treading phase is located on the left- 
hand-side of the dashed lines. The solid lines represent the 
tumbling-to-swinging transitions. 



6* < for larger 7*. Near the boundary of these two bi- 
furcations, a second range of stable fixed points appears 
between the saddle and unstable fixed points at < 0. 
However, these stable points at 6* < will likely disap- 
pear when the shape degrees of freedom in the vorticity 
direction are taken into account. 

Fig. [4] shows the full phase diagrams, both for varia- 
tions of the internal and the membrane viscosities. With 
increasing reduced volume V* the smallest shear rate 7*3, 
at which swinging can be observed, decreases, since the 
dependence of the energy F on becomes more shallow, 
compare Fig. [S] Also, the region of stability of swinging 
shifts to higher viscosity (Tyin or ?7inb) with increasing V* . 
On the other hand, for fixed V* , 7*3 is a monotonically 
increasing function of t]*^^ or 77^*^^ . This increase is found 
to be much less pronounced for rf^^ . The reason is that in 
our derivation of Eq. ([T]), we have neglect the Vj*^^ depen- 
dence. In analogy with Eq. (I7|), it is plausible to expect 
that in a more detailed calculation, a linear combination 
of and 77;*^ will appear in the both h and F; . If this 
^mb dependence were taken into account, the 7*^ lines in 
Fig. |4ja) will curve upwards more strongly; this is con- 
sistent with the larger ryj^,^ dependence observed in the 
simulations in Fig. [2] 

In the experiments of Ref. 4', the swinging (tumbling) 
motion was observed for 7* ~ 17 (4.5), = 6 (8.4), and 

- 0.1 



(calculated from r^mb ~ 10~^Ns/mof Ref. |22|) . It is men- 
tioned that the swinging motion is seen in particular in 
the close vicinity to the tank-treading-to-tumbling transi- 
tion. Thus, the experimental data confirm our predicted 
phase diagram. Fig. [U very well. Furthermore, the oscil- 
lation amplitudes also show good agreement, see Fig.O 

In summary, we have studied the oscillatory motion of 
fluid vesicles in simple shear flow. We have developed a 
simplified model for the ellipsoidal fluid vesicle, which ex- 
plains the simulation and experimental results very well. 
This model could be extended in the future to study the 
coupling of different types of oscillation mechanisms, like 
membrane shear elasticity and viscosity, in elastic cap- 
sules or red blood cells. 

Wc thank R. Finken (Stuttgart) for his help to rederive 
Eq. 11]). GG acknowledges support of this work through 
the DFG priority program "Nano- and Microfluidics" . 
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